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1 Introduction 



The Next-to-Minimal Supersymmetric Standard Model (NMSSM) [1] provides the simplest solution 
to the fi problem of the MSSM [2]. Its phenomenology can differ in various respects from the MSSM. 
Notably, as emphasized in [3], a light CP-odd scalar A\ can appear in the Higgs spectrum. 

In the case where the CP-even Higgs boson H decays dominantly into a pair of CP-odd 
scalars [3-10], LEP constraints on CP-even Higgs masses [11] are alleviated considerably [5-10]. 
For niA 1 < 10.5 GeV, where the A\ decay into BB is forbidden, this scenario could even explain 
the 2.3 a excess in the e + e~ — > Z + 2b channel for M<ib ~ 100 GeV [6] (where the two b quarks would 
result from a CP-even Higgs H with M H ~ 100 GeV and a branching ratio B (H -► bb) ~ 0.08, 
but B{H — > AiAi) ~ 0.9). Also at hadron colliders the search for CP-even scalars would be par- 
ticularly difficult [4,6,8-10] if they decay dominantly into A±Ai with itia 1 below the BB threshold. 
In this case, however, the A\ can have important effects on T decays [12-18]. Notably a Super 
B factory can then play an important and complementary role [19] via its potential sensitivity to 
T — > 7^4i decays. 

Whereas a light A\ is also possible in the MSSM with a CP-violating Higgs sector [20,21], 
scenarios with more than one gauge singlet [22], little Higgs models and non-supersymmetric two 
Higgs doublet models (see [23] for an overview) , we concentrate subsequently on the simplest version 
of the NMSSM with a scale invariant superpotential. 

T — > 7A1 decays have been investigated in the NMSSM before in [16] for m^j < 9.2 GeV, where 
the signal relies on a narrow peak in the photon spectrum. Recent results from CLEO on radiative 
T(l£) decays [24] (assuming a A\ width < 10 MeV) constrain this domain of strongly. One 
of our aims is to translate the CLEO results into constraints on the NMSSM parameter space (see 
also [25]) as Xd, the (reduced) coupling of A\ to b quarks, and to compare them with constraints 
from B physics [26-28] and the anomalous magnetic moment of the muon [25,29]. 

For niA 1 > 9 GeV, various corrections to the T — > jA\ decay rate become relatively large and 
uncertain [30], which makes it difficult to translate experimental constraints on the decay rate into 
constraints on Xd- Consequently, this region of n%A 1 is hardly constrained by CLEO results. 

Very recently, a CP-odd state with a mass of about 9389 MeV has been observed in Y(3S) 
decays by BaBar [31], showing up as a peak with a significance of 10 standard deviations in the 
photon energy spectrum. At first sight, this state can be interpreted as the long-awaited rj^lS). 
However, in the presence of a CP-odd Higgs with a mass in the same region, the observed mass has 
to be interpreted as an eigenvalue of a 2 x 2 mixing matrix, and would differ correspondingly from 
m Vb0 , the mass of the 77b (15) in the absence of a CP-odd Higgs (m^ b0 is now one of the diagonal 
entries of the mass matrix) [12]. In this case, a second peak in the photon spectrum should possibly 
be visible; however, the search for such a second peak would require a dedicated consideration of the 
various background contributions (notably from the ISR and x&j(2P)), which should be performed 
in the future. 

First, this mixing effect could explain the fact that the observed mass is somewhat lower than 
expected, if uia 1 is somewhat above m VbQ . Second, the off-diagonal element of the mass matrix 
can be estimated and turns out to be proportional to Xd [12,17]. Assuming a reasonable range for 
m 1]b0 , the observed value of ~ 9389 MeV for one of the eigenvalues implies an upper bound on Xd 
as a function of m^, which is, however, particularly strong only for mA 1 ~ 9389 MeV and will be 
derived below. 

At present, a direct detection of a CP-odd Higgs with a mass in the particularly interesting 
region 9.2 GeV < m,A 1 ^ 10.5 GeV via a peak in the photon spectrum seems to be quite difficult. 
Fortunately, an alternative signal for an A\ state below the BB threshold can be a breakdown of 
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lepton universality (LU) in T — > decays (via an intermediate A± state), since A\ would 

decay practically exclusively into t + t~ [13, 15, 17]. Note that, to this end, the photon does not 
have to be detected. Present tests of lepton universality in T(15), T(25) and T(35) decays [32] 
have error bars in the 5-10% range. Remarkably, however, a general trend (at the 1 a level) seems 
to point towards a slight excess of the r + r~ branching ratios, as expected in the presence of a A\ 
state. Another aim of the present paper is to investigate corresponding sensitivities of forthcoming 
experimental data, assuming a possible reduction of the errors to the 2% range. 

The layout of this paper is as follows: In section [2] we review the domains of the NMSSM 
parameter space which lead to a light CP-odd Higgs with strong couplings to down type quarks 
(and leptons). In section [3] we derive constraints on the NMSSM parameter space from recent 
CLEO results, using quite conservative estimates for the corrections to the T — ► ^A\ decay rate 
which lead to quite conservative upper bounds on Xd as a function of m J 4 1 . (These upper bounds 
on Xd will be included in future versions of the code NMSSMTools [33].) In section 0J we discuss 
the mixing effects of A% with r]b(nS) following [12,13,17]. In section [51 we derive constraints on 
Xd from the measured f]b obs mass by BaBar and from (conservative) assumptions on m Vb0 , and 
discuss quantitatively the possible mixing- induced shift of the measured t]b obs mass. In section [6] 
we compare these CLEO and BaBar constraints with constraints from LEP, B physics and the 
muon anomalous magnetic moment. This analysis is performed with the help of the updated 
NMSSMTools package. In section [7] we reconsider A± masses between 9.2 and 10.5 GeV and show 
that (for less conservative estimates of the corrections to the T — > jAi decay rate) a breakdown of 
lepton universality in T — » (7) l + l~ decays can become an important observable for the detection of 
a CP-odd Higgs in this mass range. We present formulas for the relevant branching ratios including 
possible A\ —rjb(nS) mixings, and study future sensitivities on Xd from lepton universality breaking. 
Section [8] contains conclusions and an outlook. 

2 A light CP-odd Higgs in the NMSSM 

In this section we show that the parameter space of the NMSSM can accomodate a light CP-odd 
Higgs, which is strongly coupled to down-quarks and leptons (see also [3,6,7]). We consider the 
simplest version of the NMSSM with a scale invariant superpotential 

W = XSH u H d + ^kS 3 + . . . (1) 

and associated soft trilinear couplings 

V soft = (\A x SH u H d + ^A K S 3 ) + h.c. + ... (2) 

in the conventions of [33]. A vev of the singlet field s = (S) generates an effective //-term, and it 
is convenient to define also an effective 5-term: 

fj, eS = As, B cS = A x + ks . (3) 

The Higgs sector of the NMSSM contains six independent parameters, which can be chosen as 

A, k, A x , A K , tan/3, fi cS . (4) 

In the NMSSM, two physical pseudoscalar states appear in the spectrum, which are superposi- 
tions of the MSSM-like state Amssm (the remaining SU(2) doublet after omitting the Goldstone 
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boson) and the singlet-like state As- In the basis {Amssm, As), the 2x2 mass square matrix for 
the CP-odd Higgs bosons has the following matrix elements [33] 

2 = M 2 2 = Xv(A x -2 KS ) 

sin Ip 

„^ 2 A 2 v 2 sin2/3 . t A . n 

M22 = o {A x + 4ks) - 3ksA k (5) 

where v 2 = l/(2y/2G F ). The masses of the CP-odd eigenstates Ay 2 are 

< 2 4[ M ii+4tAM 2 ] (6) 

with AM 2 = - M 22) 2 + 4 ( M i 2 ) 2 - 

The lighter CP-odd state Ay can be decomposed into (Amssm, As) according to 

Ay = cos O a A M ssm + sin 6 A A S , (7) 

where the mixing angle 9 a is 

M? 9 - Ml 

cosWa = ^Km^- (8) 

To a good approximation (for moderate A\ , small A K and large tan (3) , the mass of the lightest 
CP-odd Higgs boson and cos 9a can be written as [3, 7] 

~ iKS \ 2» eS B eS - 3ksA k sin 2(5 ^) ' (9) 

A«(A A -2/cs)sin2/9 

cosb'A — j -• (10) 

2fi e gB e g + 3ksA k sin 2p 

(The approximate equation for cos 9a ceases to be valid if the second term in the denominator is 
large compared to the first one.) 

The reduced coupling of the light physical Ay Higgs boson to down-type quarks and leptons 
(normalized with respect to the coupling of the CP-even Higgs boson of the Standard Model) is 
given by 

Xd = cos 9a tan f3 . (11) 

Interesting phenomena in the T-system take place for large values of Xj, i.e. large values of 
tan j3 without cos 9a being too small. (A possible enhancement of X^ [18] can occur due to the 
radiatively generated tan /3-enhanced Higgs-singlet Yukawa couplings [34]. However, in the case of 
a sizable value of cos 9a already at tree level as considered below, this effect is small.) 

At first sight eq. (fTUj) seems to imply (from sin 2(3 ~ 2/tan/J for large tan/3) that cos 9a 
decreases indeed with tan/3 - this would be the case in the PQ-symmetry-limit (k — * 0) or 
R-symmetry- limit (A K , A\ — * 0), where the second term in the denominator of (|10p tends to 
zero. On the other hand, it follows from the minimization equations of the scalar potential of the 
NMSSM (as in the MSSM), for fixed soft Higgs mass terms and n e g, that tan/3 is proportional to 
l/|^ e ff-B e ff| for large tan/3 [35], hence large values of tan/3 are associated to small values of |B e ff| 
(since |/x e ff| > 100 GeV from the lower bound on chargino masses). It is useful to replace |-B e ff | by 
the parameter 

«i=«a-w- (12) 



3 



which sets the scale for the masses of the complete SU(2) multiplet of Higgs states including a scalar, 
a pseudoscalar and a charged Higgs as in the MSSM (in our case, the corresponding pseudoscalar 
is the heavier one A2). In terms of M\, Xd can be written approximately as 



\v(A x - 2ns) 



x tan f3 , 



(13) 



M\ + 3ksA, 



and it is reasonable to examine the large tan (3 region keeping Ma fixed. 

It follows from eq. that there exist always values of A K of the same sign as A\ (typically 
both negative) where rriA 1 is small [9], while cos^ ~ 0.1 — 0.6 and hence Xd is not suppressed. This 
requires a moderate fine-tuning of A K (or Ma)', on the other hand the authors of Ref. [16] stress 
that such values for A K , which allow a light SM-like Higgs to decay into two A\ with < 2m;,, 
correspond to the smallest degree of fine-tuning in the entire parameter space of the NMSSM. For 
corresponding values of A K , the denominator of Xd in (|13p is dominated by M\. For M\ < |/cj4 k s|, 
even larger values of cos 6a — 0.6 — 1.0 are possible while m^j remains small. In this regime, the 
approximations leading to eqs. Q, (fTOj) and (fl~3|) are no longer valid, however. 

To summarize, the following conditions can be fulfilled simultaneously in the NMSSM, which 
yield possibly observable effects in T decays: 

• mA 1 < 10.5 GeV from, e.g., appropriate values of A K ; 

• a large value of Xd, if tan/3 is large while Ma in the denominator of (|13p remains moderate. 
The numerical results in section [6] confirm the analytical estimates above. 

3 Constraints from CLEO 

Recently, the CLEO collaboration presented results on Higgs searches from T(IS') decays [24]. 
21.5 • 10 6 T(15) decays had been collected and, for the T(15) — > 7 + {A\ — > t + t~) search, the 
photon energy spectrum in events with missing energy and one identified fi^ or e ± (allegedly from 
r — > evv or t — > fjLvv) had been examined. For the A\ — > /j. + /j>~ search, both muons were identified. 

No narrow peaks (of a width below ~ 10 MeV) in the photon energy spectrum are observed 
(except for T(1S) — > jJ/^Sf — > 7/i + /i~), which allows to place stringent upper limits between 10~ 4 
and 10~ 5 on the branching ratio B (T(15) -» 7(^1 -» t + t-/^ + ^-)) for m Al < 9.2 GeV [24]. 

The B (T(IS') — > 7^1) is given by the Wilczek formula [36,37] 



where a denotes the fine structure constant and Xd is given in pip. F is a correction factor, which 
includes three kinds of corrections to the leading-order Wilczek formula (the relevant formulas are 
summarized in [30]): bound state, QCD and relativistic corrections. Bound state effects have a 
quite different behaviour for a scalar or a pseudoscalar Higgs, increasing the ratio (j!4j) by ~ 20% 
in the latter case [38-40]. QCD corrections reduce the ratio p4|) by a similar amount [41,42]. 
Relativistic corrections can generate an important reduction, and were calculated in [43]. 

These relativistic corrections depend quite strongly on the b quark mass m&, and become unre- 
liable at least for Higgs masses mA 1 above 8 GeV [43] where they can generate a vanishing (or even 
negative) correction factor F. Frequently, the approximation F ~ 0.5 for all is employed in 
the literature [26,37]. However, in order to derive conservative bounds on the NMSSM parameters 



g(T(15) -> 7 Ai) 
B(T(1S)^m + M-) 
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Figure 1: F{niA 1 ) including the bound state and QCD corrections, and a naive extrapolation of the 
relativistic corrections computed for rriA 1 "C mi at larger values of m^. Black curve: assuming 
nib = 4.9 GeV (as used later), green curve: assuming nib = 5.3 GeV. 



from CLEO results, we use in this section the smaller values of F{rriA 1 ) for larger m^, which are 
obtained by a naive extrapolation of the relativistic corrections [43]. Using the quark model value 
nib = 4.9 GeV, the resulting behaviour of F{niA 1 ) (including also the bound state and QCD cor- 
rections) is shown in Fig. [IJ according to which F vanishes (and even becomes negative, in which 
case we take F = 0) for niA 1 ^ 8.8 GeV. Correspondingly, the CLEO bounds on the NMSSM 
parameters disappear for niA 1 <i 8.8 GeV. (For larger values of nib as 5.3 GeV, F would vanish 
only for niA x ~ 9.4 GeV ~ mj as also indicated in Fig. [TJ) 

Next, in order to translate the CLEO bounds into bounds on X^niA^ using eq. (fT4"|) . the 
branching ratios B (A± — > t + t~ / /i + n,~) have to be known, which depend essentially on tan/3. For 
niAj above 2m T , B {A\ — > t + t~) varies from ~ 70% for tan/? = 1.5 to ~ 95% for tan/? = 50, 
whereas B {A\ — > is always below 10% even for niA 1 below 2m T (which implies to reconsider 

the estimates of the CLEO reach in [44]). Using the code NMSSMTools [33] for the determination 
of the B (A\ — > t + t~ I ' n + and an interpolation of the CLEO bounds [24], we show our resulting 
upper limits on as a function of m J 4 1 for two extreme values of tan/3 = 1.5 and 50 in Fig. [2l 

Actually, in Ref. [24] the total decay width of A%, r^, is assumed to be below 10 MeV. 
Although we do not believe that the CLEO bounds disappear completely in the case where Ta 1 
(which increases with and niA x ) is larger than 10 MeV, we indicate in Fig. [2] also the region at 
large X^ and niA x ^ 3.5 GeV where Ta 1 exceeds 10 MeV (depending also slightly on tan/3). In 
the updated version 2.1 of the NMSSMTools package [33] these bounds are included. 
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Figure 2: Upper bounds on as a function of for two extreme values of tan/3 = 1.5 (red 
curve) and tan/3 = 50 (black curve) using results from CLEO [24]. We also indicate as dashed 
lines the region at large and m J 4 1 > 3.5 GeV where exceeds 10 MeV (same colour code for 
tan/3 = 1.5 and 50). 



M n = r 2 W 10 2 • r ( 15 ) 

n \ om„ m — riqm - - 1 ' — 1 



4 Mixing of A\ with the r)b(nS) resonances 

In the presence of a pseudoscalar Higgs boson with a mass close to one of the different 7]b(nS) 
resonances, a significant mixing between these states can occur [12]. 

The mixing between a CP-odd Higgs and a single rjb(nS) (n = 1,2 or 3) resonance can be 
described by the introduction of off-diagonal elements denoted by <5m^ in the mass matrix [12, 17] 
(here and below we neglect possible induced rjb(nS) — r]b(n'S) mixings for n ^ n') 

m\ - im Al0 T Al0 5m 2 n 

2 _ ■ 

\o (nS) lm rib0 (nS) L r) b0 (nS) 

where the subindex '0' indicates unmixed states: itia 10 and m Vb0 ^ nS ^ (Ta 10 and T Vb0 (nS)) denote the 
masses (widths) of the pseudoscalar Higgs boson and r]bo(nS) states, respectively, before mixing. 

In Ref. [17] only the mixing of the Higgs with the ^(lS) resonance was taken into account. In 
this paper, we extend the analysis by considering the possible mixing between the Higgs and any 
of the three rjbo(nS) (n = 1,2 or 3) states. Thus three mixing angles have to be defined; however, 
only the contribution from the closest r]b(nS) state to the hypothetical At mass will be assumed to 
be significant for the mixing, i.e. only one among the three mixing angles will deviate significantly 
from zero. The generally complex mixing angle a n between the pseudoscalar Higgs vlio and an 
rjbo(nS) state is given by [17] 

sin2a n = <5m^/A^ (16) 

where 

^n = [Dl + (8mlf]^ (17) 



G 



with 

D n = (m 2 Al0 ~ m 2 Vw{nS) - im Aw T Aw + im %[)(n5) r %o(n5) )/2 . (18) 

The off-diagonal element <5m 2 can be computed within the framework of a non-relativistic quark 
potential model as 

(3m 3 \ V 2 
-g^J \R Vb (nS)m X X d . (19) 

In a non-relativistic approximation to the bottomonium bound states, the radial wave functions 
at the origin can be considered as identical for vector and pseudocalar states, i.e. Rr(nS)(^) — 
Rrib(nS)(P)> an d can therefore be determined from the measured T — ► e + e~ decay widths: 



Rr(ns)(0) \' z * r[T(n5) - e + e-} x 



1 + 7T T ' 

3ir 



(20) 



Substituting recent values for the dielectron widths from [32] we obtain |-R^(is)(0)| 2 = 6.60 GeV 3 , 
l^ 6 (25)(0)| 2 = 3.02 GeV 3 and |i? %{35) (0)| 2 = 2.18 GeV 3 , leading td 

8m\ = 0.14 GeV 2 x X d , bm\ = 0.11 GeV 2 x X d , 5ml = 0.10 GeV 2 x X d . (21) 

The A\ and rjb(nS) physical (mixed) states can be written as 

Ax = cosa n vlio + sina n r] bQ (nS) , 
r] b (nS) = cos a n r] bQ (nS) - sina n Axo (22) 

assuming cos 2 a n +sin 2 a n ~ 1, i.e. neglecting the imaginary components of a n . (Here and below we 
use the notation Ax and rj^nS) for the mixed states in order to indicate their dominant components 
for small mixing angles. Clearly, for a n ~ 90°, their dominant components would be reversed.) 

The full widths T Al and T Vb ( nS ^ of the Ax and r] b (nS) physical states can be expressed in terms 
of the widths of the unmixed states according to [17] 

T Al ~ cos 2 a n r^ 10 + sin 2 a n r %0 ( n5) , 

r % ~ cos 2 a n T Vbo{nS) +sin 2 a n r Al0 . (23) 

Finally, let us recall that the mixing of the Axo with %o states should lead to mass shifts which 
can be sizable [12, 17]. These mass shifts might have spectroscopic consequences concerning the 
hyperfine 7] b (nS) — T(nS) splitting [12,17,46] whose predictions within the SM are reviewed in [47], 
and with respect to which the BaBar result [31] on the rj b (lS) — T(1S) hyperfine splitting - in the 
absence of a light CP-odd Higgs - would be somewhat large (see the next section) . 

5 Upper bounds on from the measured rji mass, and the 
mixing-induced rjb mass shift 

The observation of an r^-like state with a mass of ~ 9.389 GeV by BaBar [31], allows to obtain 
upper limits on the reduced coupling X d as a function of the lightest CP-odd Higgs mass parameter 
m Al0 , if m Al0 is near 9.39 GeV. This follows from the fact that the measured mass squared has 



3 Similar values can be obtained from a Buchmuller-Tye potential [45]. 
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now to be considered as (the real part of) the eigenvalue of the matrix M\ ()15p . corresponding 
algebraic relations and an estimate of hadronic parameters as m^iis)- 

Subsequently we denote the "77b" mass as measured by BaBar by m obs , and the state 7760(15") 
by 7750- The observed state has now to be considered as a superposition of A\q and r]bo- Then the 
following algebraic identity holds (where 5m\ is the off-diagonal element of the matrix A4 2 f)15j) ) : 



(5ml) 2 = 

where 



7 2 

1 + 



(24) 



A A = m 2 AiQ - m 2 obs , A v = m 2 w - m 2 obs (25) 

and 

7 = m Al0 T Aw - m Vb0 T Vb0 . (26) 

Note that A^ and A v must have the same sign, which follows already from properties of 
eigenvalues of real 2x2 matrices. 

Now, if we use estimates for the parameters m^ and 7, eq. (|24p allows to obtain an upper 
bound on as a function of m Aw - First, for 7 we can assume I7I < m obs x 20 MeV (from Ta 10 , 
^Vbo ~ 20 MeV) with the result that the term ~ j 2 in (|24l) is relevant only for rriA 10 very close to 
m obs . 

For (mA 10 , m %0 ) ~ m obs (but (niA 10 , m %0 ) — m bs larger than a few MeV such that the term 
~ 7 2 can be neglected), eq. (|24p can be simplified further with the result 



(5mf) 2 ~ Aml bs (m Al0 - m obs ) (m rib0 - m obs ) (27) 



and hence, from (|2ip . 



X d ~ 125 J (m Al0 - m obs ) (m Vb0 - m obs ) GeV 1 . (28) 

To proceed further, we have to estimate m Vb0 . Most of previous estimates for correspond 
actually to m VbQ — m obs > [47,48], but subsequently we allow for 

m Vb0 - m obs = -30 ... +40 MeV . (29) 

Next we have to treat the cases 7ba 10 — m obs > and rriA 10 — m obs < separately. Starting with 
rriA 10 — m obs < 0, the maximally possible value for from (|28|) is assumed for the lowest estimate 
of 

m %o' with the result 

XT X (m AlQ ) ~ 22 y/ mobs ~ m Al0 GeV" 1 / 2 . (30) 

For rriAio ~ mobs > 0, on the other hand, the maximally possible value for Xj is assumed for the 
largest estimate of m^, with the result 

JTK) ~ 25^m Al0 - m obs GeV" 1 / 2 . (31) 

These analytic expressions for X'^ MX (m Al0 ) are fairly good approximations to the numerical 
upper bounds on Xd(m Al0 ) which can be derived from ([24^1 without the approximation ([27]) . apart 
from the region where |t71a 10 ~~ m obs\ is less than about 0.5 MeV (where (|30|) and (|3ip would 
imply X™ ax (rriA 1Q ) 0). In fact, with I7I ~ m obs x 20 MeV, one obtains X™ ax (rriA 10 ) ~ 0.6 for 
\ m A 10 ~ m obs\ ^ 0.5 MeV (see Fig. [3] below) . 
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We emphasize, however, that most previous estimates for m Vb0 correspond to m Vb0 — m b s > 
[47,48] in contrast to our more conservative assumption (|29p . These estimates can still be correct 
within the present framework, if an additional Aiq state with mA 10 — m & s > exists, which mixes 
strongly with the ijbo, reducing the lower eigenvalue of the mass matrix (I15p , The induced mass 
shift ra^ w — m b s can easily be derived from eq. (f28|) : 

X 2 , x 1 GeV 2 . , 

- m obs ~ - r 32 

1.56 • 10 4 • (m Aw - m obs ) 

which, for — m b s > 0, would unwittingly be interpreted as an excess of the "observed" 
hyperfine splitting m^ns) — mobs- For instance, an induced mass shift of m Vh0 — m b s ~ 20 MeV 
would be generated by a CP-odd Higgs with mass rriA 10 and a reduced coupling satisfying 
Xd ~ 17.7 x y/mX[ — m bs GeV~ 1//2 as, e.g., X d ~ 12 for rriA 10 — 9.85 GeV. 

Note, however, that this mecanism would imply a heavier mass eigenvalue of the mixing matrix 
()15p is not too far above m b s - This favours a heavy mass eigenvalue below 10.5 GeV, which not only 
satisfies LEP constraints but could even, as mentionned in the introduction, explain an observed 
excess of events at LEP. 



6 Comparison of constraints from CLEO, BaBar, B physics and 
the muon anomalous magnetic moment 

In addition to the constraints obtained in section [3] from CLEO and in section [5] from BaBar, 
the (m An Xd) -plane is already constrained by processes from B physics [26,27] and the muon 
anomalous magnetic moment (g — 2)„ [29]. (Upper limits on Xd have been derived by OPAL [49] 
from Yukawa production of a light neutral Higgs Boson at LEP, which seem more restrictive than 
the constraints from CLEO for tua 1 ^ 9.2 GeV. We believe, however, that the r]b(nS) — A± mixing, 
which is relevant here, depends on an additional b-b-ijb form factor, where the initial 6-quark is far 
off-shell. Since this effect has not been considered in [49], we will not consider the corresponding 
limits below.) 

In the following, we will compare the different constraints in the (X^,m J 4 1 )- and (X^, Ma)- 
planes. (In this section, rriA x is the CP-odd Higgs mass parameter denoted as rriA 10 in the mass 
matrix (fT5j) . However, the difference between m^ 10 and mA 1 would hardly be visible in the Figures 
below.) 

For this purpose we have performed a scan over the NMSSM parameter space using the 
NMHDECAY program from the NMSSMTools package [33]. NMHDECAY allows to verify si- 
multaneously the phenomenological constraints from SUSY searches, Higgs searches, B physics 
and (g — 2)^. We have varied the NMSSM parameters A, k, A\, A k , /i e g, tan/3 (the latter 
between 1 and 50), as well as the SUSY breaking gaugino, squark and slepton masses and trilinear 
couplings, keeping only points where m J 4 1 < 10.5 GeV. Then we identified regions in the parameter 
space which are ruled out by the various phenomenological constraints for any choice of parame- 
ters. In particular, LEP constraints from Higgs searches require tan/3 > 1.5 in the NMSSM, while 
constraints from (g — 2)^ lead to tan (3 > 2 for m,A 1 < 10.5 GeV. 

The various curves in the (X^, m^^-plane in Fig. [3] indicate lower bounds on Xj from various 
phenomenological constraints. We found that even for very large X^ there always exist parameter 
choices such that no region is always excluded by either the constraints from B (B — » X s j) or 
B (B + — > t + v-A. However, constraints from B (B s — ► /i + ^~) and AM g , q = d,s (shown as a green 
dashed line) always exclude a funnel for m,A 1 ~ Mg q ~ 5.3 GeV, the width of which depends on 
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m A [GeV] 



Figure 3: Upper bounds on Xd versus the A\ mass for all parameters scanned over (see text). 
Indicated are constraints from B s — * fjr fj>~ and AM q , q = d, s as a green dashed line, constraints 
from {g — 2)^ as a blue dashed line, the latest bounds from CLEO on B (T — ► jtt) as a black line 
and constraints due to the measured i]b(lS) mass by Babar as a red line. 

the loop-induced b — s — A\ coupling. This coupling being proportional to Xd, the excluded region 
broadens steadily with Xd as can be observed in Fig. El leading to the exclusion of all pseudoscalars 
with masses below ~ 6 GeV for Xd > 30. However, the CLEO constraints indicated as a black line 
are much more restrictive, apart from a narrow window around m J 4 1 ~ 5.3 GeV. 

Constraints from {g — 2)^ originate from the contribution of a light pseudoscalar, which is 
enhanced by X%. For tua 1 below ~ 3 GeV, the pseudoscalar contribution has the opposite (negative) 
sign with respect to the deviation of the Standard Model prediction from the measured value of 
(g — 2)^ [29]. This results in the exclusion of very light A\ below ~ 2 GeV for large values of Xd, 
as indicated in Fig. [3]- a region which is now also covered by CLEO constraints. 

Finally, constraints due to the measured r]b(lS) mass by Babar as discussed in section [5] exclude 
a funnel around m J 4 1 ~ 9.4 GeV, which is outside the region covered by CLEO. 

In the (Xd, Ma) -plane shown in Fig. H] one sees that, as discussed qualitatively in section [21 
large values of Xd can occur only for not too large values of Ma; here this statement can be 
verified quantitatively. No region is generically excluded by the constraints from CLEO, Babar or 
(g — 2)^, since m J 4 1 varies from 1 to 10.5 GeV for each point in this plane. The upper bounds 
on Xd from B s — > and AM g , q = d,s are indicated as in Fig. [31 and lower bounds on 

from B (B — > X s j) now appear as well at small Ma- (There, the contribution to B (B — > X s ^) 
from a charged Higgs with a mass ~ Ma has to be compensated by a contribution ~ X^ involving 
charginos or neutralinos, which requires a sufficiently large value for Xd-) 

The most important results of this section are contained in Fig. [31 which shows that the com- 
bined present constraints rule out most of the region where mA 1 ^ 8.5 GeV - except if Xd is 
sufficiently small - whereas mA 1 ~ 8.5 — 10.5 GeV remains an interesting region in parameter space 
allowing for large Xd- The corresponding necessary values of Ma can be deduced from Fig. |H 
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Figure 4: Bounds on Xd versus Ma for A% masses below 10.5 GeV. Indicated are upper bounds 
on Xd from B s — > /U + /U~ and AM,, q = d, s as a green dashed line, and lower bounds on X^ (for 
m A < 200 GeV) from B {B -> X s j) as a blue line. 



7 Possible lepton universality breaking 

As pointed out in [13], one manifestion of the existence of a light CP-odd Higgs boson of mass 
around ~ 10 GeV could be a breakdown of LU (lepton universality) in T decays, if the (not 
necessarily soft) radiated photon escapes undetected in the experiment, or simply is not specifically 
searched for in the analysis of events. (The leptonic width is, in fact, an inclusive quantity with a 
sum over an infinite number of photons.) Higgs-mediated T decays would lead to an excess of its 
tauonic branching ratio (BR), which can be assessed through the ratio 

„ B TT — Bff B TT , N 

n T/l = TT =ir~ l ( 33 ) 

ott oil 

where B TT denotes the tauonic, and Bu. the electronic [t = e) or muonic (£ = fi) branching ratios of 
the T resonance, respectively. A statistically significant non-zero value of 1Z T /£ would be a strong 
argument in favour of a pseudoscalar Higgs boson mediating the process. 

In Table Q] we summarize the current situation of LU obtained from [32]. As already mentionned 
in the introduction, a ~ 1 a effect seems visible in most cases, leading to an overall (positive) ~ 2 a 
effect. 

Subsequently we intend to estimate the possible amount of LU breaking in the NMSSM with 
mA 1 in the 9 - 10.5 GeV range and large X^ in order to verify whether it can be assessed experi- 
mentally [17], i.e. whether it can be of the order of the few percent. 

In principle, also pure SM channels may yield an apparent breaking of LU in T decays. In 
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B(e + e~) 


b(a* + m-) 


B{t+t-) 


R T/e (nS) 




T(15) 


2.38 ±0.11 


2.48 ±0.05 


2.60 ±0.10 


0.09 ± 0.06 


0.05 ± 0.04 


T(25) 


1.91 ±0.16 


1.93 ±0.17 


2.00 ±0.21 


0.05 ±0.14 


0.04 ±0.06 


T(35) 


2.18 ±0.21 


2.18 ±0.21 


2.29 ±0.30 


0.05 ±0.16 


0.05 ±0.16 



Table 1: Measured leptonic branching ratios B(Y(nS) — > l€) (in %) and error bars (summed in 
quadrature) of T(1S), T(25), and T(3S) resonances. 

fact, all decays through an intermediate pseudoscalar (e.g. rjb) state should break LU due to the 
leptonic mass dependence of the amplitude as a consequence of helicity conservation. However, 
such processes (mediated by a two-photon loop or a Z°-boson in the SM) provide contributions to 
the branching ratios well below 1% and can be dropped. On the other hand, phase space should 
suppress the tauonic BR (by about 0.5%) with respect to the electronic and muonic modes. 

In the absence of mixing between the Higgs boson A\ and the % resonances, the relevant 
contribution to TZ T /£ would originate exclusively from the first diagram a) in Fig. Assuming a 
BR(Ai — » t + t~) of 100% (see, however, below), R T /£ would be given by the Wilczek formula (j 14|) : 




b 



Figure 5: Process e + e~ — > T — ► 7 t + t~ with a) pseudoscalar Higgs, b) % (after mixing) as 
intermediate states; c) mixing diagram. Diagrams similar to a) and b) could be drawn for the 
two-gluon decay mode yielding hadrons in the final state instead of taus. 

We recall that here we are interested in A\ masses above 9 GeV, in which case the very con- 
servative (small) estimate in section [3] of the correction factor F in (|34p ceases to make sense. The 
aim of section [3] was to derive conservative upper bounds on X^; here, however, we aim at realistic 
estimates of 1Z T /e- To this end, the tree-level expression (|34"|) corrected by a constant factor F = 1/2 
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should provide an acceptable approximation [26,37]. Then, the CLEO bounds exclude large values 
of Xd for m Al up to roughly 9.2 GeV. 

In the presence of mixing between the Higgs boson A± and one of the the ijb resonances as in 
section HI both eigenstates defined in eqs. (|22|) would contribute to TZ T /e- In order to give the T 
branching ratios into these mixed states, it is convenient to define the ratio 

5o(n ' n)= g(T(nS)->„V) ' (35) 

where the branching ratio for a Ml transition between T(nS) and r/bo( n ') states (n' < n) is given 
by [48] 

B (T(„5)^™„(n'5))=^(^) 2 %^^ (36) 

k is the photon energy (depending on the mass difference mx(nS) ~~ m r] b (riS))i ^n'n denotes the 
final and initial wave functions overlap, I n / n = {f n '\jo(kr /2)\i n ) , where jo is a spherical Bessel 
function. I n / n is numerically close to unity for favoured transitions (n = n') but much smaller for 
hindered (n ^ n') transitions. As stressed in [48], however, the considerably larger photon energy k 
in the latter case could compensate this reduction, leading to competitive transition probabilities. 
Below we set I%2 = 0.057 [48] and I13 = 0.017. (The latter value is required in order to reproduce 
the experimental value B (T(35") — > 7 hadrons) = 4.8 x 10 -4 found by BaBar [31].) 

In terms of Rq and So defined above, the T branching ratios into the mixed states A\ and 
rft,(n'S) are given by (neglecting interference terms, and normalized w.r.t. the branching ratios into 
l + l~ = or e + e~) 

B{T(nS)^ 1 A 1 ) , D , „, 2 



B(T(nS) -> 1+1-) 
B (T(nS) -> jrjbjn'S)) 
B(T(nS) -» 1+1-) 



cos a k Rq + sin a k So(n,k) 
sin 2 a n i Rq + cos 2 a n > So(n,n') (37) 



where, as stated above, we assume that at most one possible mixing angle a k is nonvanishing. (For 
a given value for mA 10 , the index k is given by the state rjbo(kS) whose mass is closest to m Al0 - The 
T decays into the remaining unmixed r)t,(n'S) states with n' 7^ k are still described by eq. (j36f) . ) 

Next we assume that the mixed states A\ and rj^nS) decay into t + t~ only via their Aiq 
component. Then we obtain, using eq. (|23j) for the full widths of the mixed states, 

cos 2 a k T A 



B{A 1 ^t + t-) = B{A 10 ^t + t-)x— . 

cos^ a k T Aw + sm z a k T Vb0 ^ 

B( Vb (nS)^T+r-) = B(A 10 ^t+t-)x - sin2a "^ ? (38) 

cos^ a n F Vb0 (nS) + sin a n T Al0 

i.e. B{r]b{nS) — > t + t~) vanishes for n 7^ k. Finally we obtain for lZ T /g (for a given T(nS), and 
assuming B (A 10 -► r+r") = 90%) 

= ^(T(n5)^ 7 ^) B(TK). 7 #5)) 

- B(T(n5)^Z+/-) l l j+ B(T(nS)^l+l-) * B^k!,) ^ t r ) 

cos 2 a fc T^io 



0.9 (cos 2 a k Rq + sin 2 a k So(n, fc)) 



cos 2 ctfc r^io + sin a fc T 



Vbo( k S) 
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+ 0.9 (sin 2 Qfc Rq + cos 2 a k So(n, k)) 



sin 2 a k T Al 



(39) 



cos 2 a k r %o(fcS) + sin 2 a k T Aw 

), or a k is given by the (supposedly 



where either a k = (if none of the states ^(kS) mixes with Aiq ) 
only) mixing angle, whose choice and value depend on m^ 10 . 

In the set of plots of Fig. [6J TZ T /£ is shown for T(15), T(2S) and T(35), respectively, as a 
function of rriA 10 , using the formulas of section 4 for the determination of the relevant mixing 
angle. We assume tentativelj@ r^ M(1S) . = T VbQ , 2S ^ = r^ 60(3S) = 5 MeV and Xd = 12. (As argued in 
section 5, values for m Al ~ 9.4 ± 0.2 GeV are actually ruled out for this value of Xd- However, 
for lower values of Xd - leading to correspondingly lower values for lZ T /£ - the forbidden window 
for m Al becomes smaller. Moreover, phenomenologically interesting values for m Al <; 9.6 GeV as 
discussed at the end of section 5, are seen to generate interesting values for TZ T /e-) 



a)Y(lS) 
r[^ b0 ]=5 MeV 



0.16 
0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 






b) Y(2S) 




r[7 7b0 ]=5 MeV 











9.4 9.6 9.8 10 

m A ,„(GeV) 



9.2 9.4 9.6 9.8 10 10.2 10.4 

m Al0 (GeV) 



Figure 6: 1Z T /i versus the pseudoscalar Higgs mass for a) T(IS'), b) T(25), and c) T(35) decays 
using Xd = 12, m Vb0 (is) = 9.389 GeV [31], assuming m Vb0 (2S,3S) = 9.997,10.32 GeV respectively, 
and ^ Vb0 (is,2S,3S) = 5 MeV. The contributions from H^j e are indicated as dashed green lines, the 
contributions from 7^?% as dotted black lines, and their sum lZ T /£ as solid red lines. Larger (smaller) 
values of Xd obviously yield higher (lower) values for R T /(. 

The contributions from (dotted black line) yield the expected bumps around the respective 
7]b mass values, where the mixing angle becomes large. Conversely, the contributions from 
(dashed green line) show dips at both m,A 1 = m Vb (is) = 9.389 GeV and mA 1 = m Vb (2S) — 10 GeV, 
since they become reduced by the mixing. (The expected peak or dip at rriA 1 = m Vb (2S) — 10-3 GeV 
in Fig. [6b) is in fact invisibly small.) 

4 The following expression (asymptotically valid for very heavy quark masses) can be used to estimate the rj b (nS) 
full width from the experimentally known r/ c (nS) full width: V Vb (nS)/F Vc (nS) — (m b /mc)[as(2m b )/a s (2'm c )] 5 [50]. 
The fifth power of the a s ratio yields however a large uncertainty to the prediction. On the other hand, theoretical 
predicitions based on the expected ratio of the two-photon and two-gluon widths range from 4 to 20 MeV [51] 
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Figure 7: R T /i versus the pseudoscalar Higgs mass for T(3S) decays using Xj, = 12 and cl) 
T Vb0 = 10 MeV, c2) T Vb0 = 15 MeV, respectively. 

The higher values of lZ T /£ and the higher reach in m,A w obtained for the T(2S) and T(3S) 
(due to the dominant Wilczek mechanism of Fig. [5h.)) allow us to conclude that radiative decays 
of the latter resonances look more promising than the T(15) decays, allowing for the experimental 
observation of LU breaking (at the few percent level) at a B factory. This result is important for 
future tests of LU [19]. 

In order to study the effect of our assumption on the width T Vb0 , we present in the set of Fig. [7] 
TZ T /£ for the T(35) resonance setting r^ b0 =10 and 15 MeV, respectively. One can observe a slight 
overall decrease of TZ T /e for larger T Vb0 , as expected from eqs. (j38j) and (f39l) . 

Concerning future measurements, we assume tentatively that a combined statistical and sys- 
tematic error (summed in quadrature) of 2% for R T /i is achievable at a (Super) B factory. Fig. [8] 
shows the foreseen 2 a (95% CL) limits (green region) for testing LU using T(3S) decays for a A\ 
mass ranging in the interval 9-10.3 GeV. 

An observation of lepton universality breaking in T decays should lead to a careful search for the 
quasi-monochromatic photons shown in Figs. [5] in a sample of events firstly selected and enriched 
using 1-prong tauonic decays and requiring missing energy (neutrinos). Let us recall that there 
could be two nearby peaks corresponding to two physical eigenstates. However, large A\ or % 
widths might invalidate this search method: if the itia 1 and m Vb masses were not too different (i.e. 
less than 50 MeV), the two peaks might not be resolved experimentally but yield a broader peak 
than expected. This conventional search has been unsuccessful so far in T(IS') decays, but can 
(should) be extended to the yet unexplored radiative decays of the T(25*, 35") resonance into r's, 
according to the proposal in [15, 17]. 

Finally, a light CP-odd Higgs which mixes with one of the % states would also decay hadronically 
and eventually be visible by requiring four or more charged tracks together with a photon in 
radiative T(2S, 3S) decays, a criterium used by Babar [31] for their discovery of the r]b(lS). Whereas 
the Babar data can possibly be used to put bounds on such an additional state for certain ranges 
of its mass, we believe that a serious search for such an additional state in the inclusive photon 
spectrum from radiative T decays would require a detailed treatment of the various background 
contributions, which is beyond the scope of the present paper. Corresponding investigations are 
clearly another interesting task in the future. 
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8 Conclusions and outlook 

In this paper we have summarized constraints from and perspectives of various processes related 
to a CP-odd Higgs boson with a mass rriA 1 below the BB threshold of 10.5 GeV. Apart from m^, 
these phenomena depend essentially on its reduced coupling to 6-quarks. Within the parameter 
space of the NMSSM, relatively large values of Xd are possible (for sufficiently large values of tan (3). 

We have compared present constraints on the mA 1 — X& plane from B meson physics, the 
anomalous magnetic moment of the muon, LEP, and recent results from CLEO and Babar. In 
spite of the conservative approach towards the bound state corrections to the Wilczek formula the 
most stringent constraints originate - not astonishingly - from the dedicated (negative) searches 
by CLEO for < 8.8 GeV allowing, however, for substantial values of X& provided that 
8.8 GeV < m Al < 10.5 GeV 

Given the possible explanation of the 2.3 a excess in searches for a CP-even Higgs boson at 
LEP [6], this allowed mass range is of particular interest. We emphasize again that the interval 
9.4 GeV < niA 1 < 10.5 GeV can also have an effect on the rjb(lS) mass as measured by Babar via 
mixing, and explain the possibly excessive T(1S) — rj^lS) hyperfine splitting. 

Such a scenario can and should be tested at presently running B factories, and/or a future 
Super B factory. Obvious search strategies consist - as already performed - in radiative T(nS) 
decays into both tauonic and hadronic final states, keeping an eye on possible close (but separate) 
peaks in the photon spectrum. In addition, violation of lepton universality in inclusive radiative 
T decays can be a signal for an additional CP-odd Higgs. We have clarified that corresponding 
visible signals are well within the reach of future precision experiments. 

These searches are complementary to Higgs boson searches at colliders (like the LHC) where it 
is quite doubtful at present whether a light CP-odd Higgs decaying dominantly into t + t~ could 
be seen. On the contrary, such a CP-odd Higgs can render searches for the lightest CP-even Higgs 
boson h very difficult, if it decays dominantly into h — > j4i^4i — > t + t~t + t~ . Hence, searches at B 
factories are possibly our only windows into the light Higgs sector, if such a scenario is realized. 
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